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p^ ■ Abstract 

1/-^ I A new and simple procedure to measure the angle 7 from B^ — > Tr^vr+vr" and 

^D ■ B — > K ir^TT^ decays using SU(3) symmetry is presented. It is based on a full 

f^^ . Dalitz plot analysis of these decays. All diagrams, including strong and electroweak 

^^ I penguins, are considered in the procedure. The method is also free from final 

^^ ' state interaction problems. The theoretical error in the extraction of 7 within the 

-«.^ ■ method should be of the order of 10° or even less. Taking into account the B- 

iu: 
> 

X 



meson production in the first generation of B factories and recent measurements 
from CLEO, this method could bring the best measurement of 7 in the next years. 
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For the next years, large accelerator facilities known as B factories are going to be 
running. The goal is to produce a large amount of B mesons because their decay should be 
sensitive to CP violation |l], Q. As a consequence, one hopes to measure the three Cabibbo- 
Kobayashi-Maskawa (CKM) angles a = arg{-VtdVtl/VudV:^), (3 = arg{-VcdV;jVtdVtl) 
and 7 = arg{—VudV*f,/VcdVX) to check the standard model predictions. It is not yet clear 
which precise measurements will allow a clean extraction of a and 7. Concerning the 
last angle, many interesting methods have been proposed so far but all of them suffer 
from either experimental or theoretical problems. On the one side, theoretically clean 
procedures based in B ^ DK decays § suffer from important experimental difficulties 
and would demand about 10 years of data taking in order to extract 7 with an error, at 
least, of the order of 15°[|ll. On the other side, other methods based in B -^ tttt, B — > 71K 
01 B ^ KK decays 0] have theoretical uncertainties which would imply systematic errors 
in the extraction of 7 as large as 20° @. It is then worth looking for new experimental 
procedures to extract 7 with smaller uncertainties. 

Methods proposed so far to measure 7 are based in the study of branching ratio 
asymmetries in two body decays. In a previous letter we showed that three body decays 
could be a more interesting tool to extract CP violating phases]^. The idea is to make 
use of the fact that Dalitz plot analysis of a three body decay gives a direct measurement 
of amplitudes instead of branching ratios. In other words, one has a direct access to the 
phase of a given process. In ref. |^ we have illustrated this general remark applying it to 
the extraction of 7 studying the decay B^ -^ tt^^tt+tt^ . We showed that, due to this direct 
experimental access to phases, 7 could be extracted with a smaller statistical error than 
the usual methods based in two body decays existing in the literature. Unfortunately, the 
method suffers from one difficulty also existing in two body decays, i.e., penguin pollution: 
the bigger is the unmeasured penguin contribution compared to the tree one, the larger 
is the systematic error of the method. 

In this letter, making a full use of Dalitz plot analysis features, we present a first pro- 
cedure to extract 7 using three body decays which takes penguin contributions explicitly 
into account. The method is based in a combined study of two pairs of CP conjugated 
decays B^ —>■ ir^n'^n' and B^ — > K^n^-n-^ related by SU(3) symmetry. It is free from 
final state interaction and electroweak penguin problems. The theoretical error of the 
method is due to SU(3) breaking and the uncertainty about charm penguin diagrams. 
Considering present expectation values for 7 and f3, we estimate this error to be of the 
order of 10° or even less. This value is smaller than that of other methods experimentally 
accessible within a few years. 

Let us start the description of the method by presenting the main features of Dalitz 
plot analysis in our context. Three body decays of heavy mesons seem to be dominated 
by intermediate resonant decays[§]. All these partial intermediate channels, together with 
the direct non-resonant channel interfere to give the same — detected — three body final 



state. The Dalitz plot analysis is a powerful tool that yields a clear separation of all these 
intermediate channels; moreover, it brings a direct measurement of the amplitude, i.e., 
magnitude and phase, of the contribution of all the intermediate processes. As phases are 
always measured with respect to a given one, in order to extract a weak phase with this 
method, one needs at least two distinct channels with different weak phase 0. 

This is the case for the two pairs of decays B"^ -^ tt'^tt+tt^ and B"^ —>■ K^tt^ti^ con- 
sidered here, thanks to the presence of Xco as a possible intermediate resonance. Indeed, 
the contribution A^B"^ — ^ Xco'P'^) x -^{Xco -^ tt+tt") (where V = n or K) has no weak 
phase while for the other possible channels — as for example those mediated by a p° or 
/o resonances — the tree contribution has weak phase 7. 

Final state interaction problems do not affect this method. One can easily get rid of 
them by choosing a contribution mediated by an isospin resonance p, ^ — such as an 
/o or a possible a resonance. Indeed, a B^ — *> foT"^ decay for example, proceeds through 
a unique isospin amplitude. Thus, the method presented in this letter is completely free 
from final state interaction difficulties. In the following, we will present the method using 
aB^ ^ foV^ decay. 

Besides tree contributions driven by the weak phase 7, A{B^ -^ fo'P^) has also 
penguin contributions with different weak phases. For V = tt, top penguin diagrams have 
weak phase (3 and charm penguin ones have no weak phase; on the other side, for V = K, 
both top and charm penguins have no weak phase. Thus, for the four decays, the /o 
amplitudes are written as 

A2 = Te^(^^-^) + Pte'^^Pt+l^^ + P^e'^p^ (2) 

A3 = T'e'^^T'+i) + p'eiSp' (3) 

A^ = T'e'^^T'-/) ^P'e^^p' (4) 

where A^ = A{B+ -^ /o7r+), A2 = A{B- -^ /ovr-), A3 = A{B+ -^ UK+) and 
A4 = A{B^ -^ j'qK^). In the expressions above, 5t,5t',5p^,5p^ and 5 pi are strong 
(CP conserving) phases, T and T' contain both tree and color suppressed contributions, 
Pt includes all strong and electroweak penguin diagrams with weak phase /3, Pc includes 
all those with no weak phase and P' includes all penguins. In other words, Eqs. (|l|-§) 
include all kind of diagrams contributing to these decaysfl^. Some of these decays can 
be related by SU(3) symmetry, as we will see below. 

The left-hand sides in Eqs. (|l|^; i-e. Aj, {i = 1, ...,4), are directly measured complex 
numbers, that is, 8 real independent quantities. This is one of the main claims of this 
letter: Dalitz plot analysis brings more independent measurements than usual two body 
branching ratios; one thus has more information available for each decay that can be used 
to treat penguin and other pollutions. 



As a first step, we will exclude from our analysis charm penguin contributions which 
appear only in the first two equations. Note that this is not a strong assumption. Indeed, 
in B^ -^ 'jt'^tt'^ti' decays, as in i? — > vrvr, one expects penguin contributions to be of the 
order of 20% of tree ones0]. Moreover, charm penguin amplitudes are expected to be not 
larger than half of top ones|p.l||. Thus, in order to resolve in 7 in Eqs. (|1|-§1), the influence 



of charm penguin is expected to be small. In spite of this, in the last part of this letter, 
we will explicitly study their influence in our results. 

Now, we make use of SU(3) flavor symmetry. All contributions included in T' are 
identical to those in T except for an s quark replacing a light one. (Note that one cannot 
make the same assumption for the penguin sector because in Pt there are two different 
topologies while in P' there is only one.) One can then write, 

T' = XT (5) 

S't = St (6) 

where A = Ks/Kd ~ 0.2. 

The validity of these assumptions will be discussed later in this letter. 

One then has a system with 8 real equations and 8 unknown quantities. After a 
simple algebra to eliminate unwanted quantities, one can reduce the set of 8 equations to 
a simpler system containing the desired variable 7: 

arg{A3 - A4) = arg{e'^^^+^^ - e^^"^^"^)). (8) 

This system, containing 3 real equations, allows us to obtain 7, (3 and 6t- At this 
stage, one could expect this method to provide not only a measurement of 7 but also an 
independent measurement of (3. Unfortunately, this is only true if Eqs. (J^) and (||) were 
exact. This is not the case due to the theoretical assumptions made above. The question 
is how does any error in these equations propagate to the actual extraction of 7 and /?. 

We found that /3 is too sensitive to small uncertainties in these equations. This is 
simply due to the fact that /3 is only present in in Eqs. (|I|) and (|]), in a term which is 
small. Thus, small uncertainties in coefficients Ai 2 imply large uncertainties in (3. 

As a consequence, the method is not well suited to providing an independent measure- 
ment of (3. Fortunately, in the next years, the phase (3 will be known within an error of a 
few degrees. Our strategy is thus to consider /3 as a known parameter. We then have two 
independent variables, 7 and 6t with only two equations, e.g., real and imaginary parts 
of Eq. (0) ^ . 



Let us now study the theoretical errors of the method and their influence in the 
extraction of 7. The sources of systematic errors are 1) the validity of SU(3) symmetry 
assumptions and 2) charm penguins. 

SU(3) symmetry is only assumed for those terms carrying the weak phase 7, i.e., 
diagrams included in T. This point represents another important issue of this method: 
one does not need to make SU(3) assumptions about penguin contributions. The two main 
contributions to T are tree and color suppressed diagrams. Factorization corrections to 
exact SU(3) symmetry in both diagrams transform Eq. (^ in 

r = ^XT (9) 

In 

where fx and /tt are the kaon and pion decay constants, respectively. As fx/ fn ~ 1-2, 
this represents a 20% correction with respect to the assumption of exact SU(3) symmetry 
of Eq. (|^). The importance of non-factorizable corrections to Eq. (|^) is not yet clear 
but there are usually assumed to be small. Both theoretical [0] and experimental fl^ 
educated studies seem to favor non-factorizable corrections not larger than 10%. 

We have studied numerically the influence of this uncertainty when solving Eqs. (|^. 
For this, after changing (j^) to @ we have assumed an extra 10% theoretical uncertainty in 
the relation between T' and T and studied the propagation of this error in the extraction 
of 7. More precisely, we have first assumed a given set of values for the various parameters 
entering the right-hand-sides of Eqs. (0-§) — including 7. With them, we have calculated 
the quantities Ai, {i = 1, ...,4), using Eqs. (|l|-§) considering explicitly a 10% correction 
to (y). Then we have solved the system of Eqs. (|^ — which assumes that there is no 
such a correction — to find 7. We have finally compared the latter with the originally 
assumed value for 7. 

We found that the amount of this error in extracting 7 only depends on the actual 
value of 7 and f3, as shown in Fig. 1. This error was found to be independent of the 
actual values of T, Pt, P', 6t, Sp^ and 6p'. Fig. 1 shows that the smaller j3, the larger the 
errors. When (3^0, the error diverges. This is simply because, when /3 is exactly zero, 
Eqs. (|l|-§) with Pc = 0, are no longer independent; thus Eqs. (|^ only admit solutions 
when SU(3) symmetry assumption is exact. 

The other assumption related to the SU(3) symmetry, i.e., Eq. (^, has no influence 
in the result. To find this, we proceeded in the same way as described above — this time 
assuming 5^ 7^ 6t when calculating Ai. We found that the error in finding 7 due to this 
assumption is always negligible, independently of the values of T, Pt, P', 6t, Sp^, Sp', (3 and 

7- 

Let us now consider the error due to charm penguins. They would only affect Eqs. (|I|) 
and (0) since in the B^ -^ K^ti^ti^ decays charm penguins have the same weak phase 
as top ones and are thus already included in P'. This is a crucial remark, because in the 
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Figure 1: The error in the extraction of 7 as a function of 7. Sohd curves are the error 
due to the uncertainty in the SU(3) assumption, for 4 different values of /5. The dashed 
curve is the error due to charm penguins. 
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decay Pt/T is expected to be small|0; thus, as 7 is present in the T-term, 
any error in the penguin sector is less important. The actual importance of these neglected 
terms is in fact unknown at present, but a crude estimate gives|n] 0.2 < Pc/Pt < 0.5. 

The theoretical uncertainty in the value of this ratio propagates to the actual extraction 
of 7. To estimate the amount of this error, we made the same kind of numerical analysis 
as described above to the study of SU(3) symmetry assumptions. In other words, we 
have calculated the quantities Ai assuming that Pc is not zero, and then solve Eqs. (|^) 
— where Pc is not present. We found that in this case, the error in the extraction of 
7 does not depend on the actual value of /5, but it does depend on the value of 7, as 
shown in the dashed curve of Fig. 1. As in the SU(3) case, these results are independent 
of the actual values of 6t, Sp, 6p and P'. Nevertheless, the dashed curve in Fig. 1 does 
depend on the actual value of Pt/T. More precisely, the error is proportional to the value 
of this ratio. Results presented in Fig. 1 have been obtained using the expected values 
of Pt/T 



0.2 |T§, PjPt = 0.35 m 



and Sp^ = Sp^ [^. These results do not change if 
the last assumption is released, and change little if Pc/Pt takes another value within the 
range 0.2 - 0.5. 

Finally, we have studied how an uncertainty in (3 propagates to the extraction of 7. 
Assuming (3 is known with an error of 3° - which is the expected error after 4 years of 



running of BaBar and Belle - this contributes with less than 1° in the error to extract 7. 
Let us now analyse the total theoretical error of the method. Adding in quadrature 
all the errors we have discussed above, we conclude the following. If (3 is 50° or more, the 
total error of the method is always small: 5° or less. If (3 turns out to have the central 
value estimated by Standard Model predictions|T^, i.e. 25° (which is also in reasonable 



agreement with recent reports from Belle and BaBar [|T^), then the error may be larger: 
if 7 is in the range 50° - 80°, as dictated by an overall analysis of the unitary trianglejlTl, 



then the error is of the order of 11°; if 7 is in the range 100° - 130°, as suggested by 
experimental constraints based in non-leptonic B decays[|l9|, then the error in extracting 
7 is of the order of 5° to 8°. 

All these values for the theoretical error are smaller that what one expects from other 
methods to extract 7, proposed so far. This is valid as far as the value of /3 is larger than 
10°. 

It is important to stress that the method is also self contained for the determination 
of its systematic error. Indeed, once the actual values of 7 and Pt/T are determined from 
the experiment, the systematic error can easily be read from the numerical study made 
above. 

Our last comment deals with the experimental feasibility of this procedure. The 
channels involved in the analysis only require the measurement of charged pions and 
kaons and do not require the necessity of tagging (since the method only deals with 
charged B meson decays). In ref. ||^ we made a simulation of B^ — > tt^tt+tt^ decays; 
we showed that the extraction of partial amplitudes - i.e., what we call Ai, {i = 1, ...,4) 
in this letter - could be done, with relatively small statistical errors, with about 1000 
reconstructed events in each Dalitz plot. Recent CLEO datapO] found that B^ -^ V^p^ 
(with V = -ft', tt) have branching ratios of the order of 10~^. Thus, assuming a total 
branching ratio of the decays B^ -^ tt^tt+tt^ and B^ -^ K^tt^tt^ of the order of 2 x 10~^ 
and a 60% reconstruction efficiency ^, this method may provide a good measurement of 
7 after about 4 years of running of the ffist generation of B machines. 

In summary, we have presented in this letter a method to measure the weak CKM angle 
7 which deals with all strong and electroweak penguin amplitudes. It exploits a general 
procedure to extract weak angles using Dalitz plot analysis in three body B meson decays. 
Here, we discuss a combined study of two pairs of CP conjugated decays, B^ — > tt^tt+tt" 
and B^ -^ K^tt^tc^ which allows the treatment of penguin amplitudes. Moreover, tree 
and penguin magnitudes and strong phases are obtained from the experimental procedure. 
The method brings a measurement of 7 free from final state interaction problems. 

The procedure has two theoretical sources of uncertainty. First, it is based in SU(3) 
symmetry assumptions, even though only for the tree sector, not in the penguin one. 
Second, charm penguins can only be included in the method by considering a model to 
estimate Pc/Pt- According to the present knowledge, charm penguins would introduce a 



small systematic uncertainty. Combining both sources of error, the total theoretical error 
of the method presented in this letter is not larger than 5° if (3 is larger than 50°; if (5 has 
the expected value of 25°, then the error is of the order of 5° - 8° if 7 turns out to be in 
the second quadrant and of the order of 11° if it is in the first one. Thus, one expects the 
theoretical error of this method to be smaller than that of other procedures proposed so 
far to extract 7. 

The method presented here should be experimentally feasible at BaBar and KEK. It 
could then provide the best knowledge of 7 in the next years. In any case, this method is 
based on a new idea and a different experimental procedure, if compared to usual methods 
based in two body decays. It will then bring an independent measurement of 7. 

References 

[1] P.F. Harrison, ed. and H.R. Quinn, ed.. The BaBar Physics Book (SLAC-R-504, 
1998). 

[2] For a general review see: Y. Nir and H.R. Quinn, in B decays, edited by S. Stone 
(World Scientific, Singapore, 1994), p. 520; M. Gronau, Nucl. Phys. (Proc. Suppl.) 
B38, 136 (1995). 

[3] M. Gronau and D. Wyler, Phys. Lett. B 265, 172 (1991); I. Dunietz, Phys. Lett. 
B 270, 75 (1991); D. Atwood, I. Dunietz, and A. Sony, Phys. Rev. Lett. 78, 3257 
(1997). 



[4] For a review see for example R. Fleischer, [hep-ph/9904313| and references therein. 



[5] M. Gronau, |hep-ph/00013171 ; and references therein. 



[6] I. Bediaga, R.E. Blanco, C. Gobel, and R. Mendez-Galain, Phys. Rev. Lett. 81, 
4067 (1998). 

[7] The use of Dalitz plot analysis to measure CP violating angles has first been intro- 
duced, in a quite different context, in: G. Burdman and J.F. Donoghue, Phys. Rev. 
D 45 187 (1992) and A.E. Snyder and H.R. Quinn, Phys. Rev. D 48, 2139 (1993). 

[8] Particle Data Group, C. Caso et a/., Eur. Phys. J. C3, 1 (1998). 

[9] J.M. Gerard and J. Weyers, Eur. Phys. J. C7, 1 (1999). 

[10] Only penguins with a quark up in the loop are excluded from the analysis of the 
B^ -^ vr'^TT+TT^ decay. Nevertheless, they are expected to be even smaller than 
charm penguins. 

8 



[11] A.J. Buras and R. Fleischer, Phys. Lett. B 341, 379 (1995). 

[12] Note that, in fact, Eq. (H) is independent of 7: its right-hand side is 
arg{2i e**^^ sin 7) = 5t + '7r/2, when 7 varies between 0° and 180°. 

[13] M. Gronau and D. Pirjol, Phys. Rev. D 62, 077301 (2000). 

[14] M. Gronau, O.F. Hernandez, D. London, and J.L. Rosner, Phys. Rev. D 52, 6374 
(1995); M. Athanas et ah, CLEO Collaboration, Phys. Rev. Lett. 80, 5493 (1998). 

[15] In fact, present estimates correspond to the ratio penguin-to-tree diagrams while 
the T-term also has color suppressed diagrams. Nevertheless, as tree contributions 
are expected to dominate, this ratio should not differ significantly. 



[16] A. Dighe, [hep-ph/9811506| and references therein. 



[17] A. Ah and D. London, Eur. Phys. C9 687 (1999). 



[18] B. Aubert et al, BaBar Collaboration, |hep-ex/0008048| ; H. Aihara et al, Belle 
Collaboration, |hep-ex/001000^ . 



[19] Y. Kwon et al., CLEO Collaboration, |hep-ex/9908039 



[20] C.P. Jessop et al, CLEO Collaboration, Phys. Rev. Lett. 85, 2881 (2000). 



